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Abstract. An analytic approach for controlling quantum states, which was originally applied to fully random matrix systems 
[T. Takami and H. Fujisaki, Phys. Rev. E 75, 036219 (2007)], is extended to deal with more realistic quantum systems with a 
banded random matrix (BRM). The validity of the new analytic field is confirmed by directly solving the Schrodinger equation 
with a BRM interaction. We find a threshold of the width of the BRM for the quantum control to be successful. 
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INTRODUCTION 

Theoretical and experimental studies of controlling quantum states have been attracting much attention because of the 
theoretical progress in the field of quantum computing [ 1 ] and of the technical developments in manipulating atomic 
and molecular systems. Various approaches have been applied to such quantum systems and its target is often "simple" 
such as one-state-to-one-state control via adiabatic passage [2]. When we consider to design quantum devices with a 
large number of states interacting with a complex environment or with short-time laser pulses, however, a multi-state- 
to-multi-state control problem is not an exception but a rule. Optimal control theory [2] and genetic algorithms |3|,|4|] 
are most successful methods to solve this kind of complicated problem, but its implementation and interpretation can 
be still difficult. 

From theoretical points of view, a "complex" quantum system may be modeled by fully random matrix systems 
with generic properties J5t] or by quantum chaos systems where their classical limit shows chaotic properties [6]. Gong 
and Brumer applied the coherent control method 13,181] to a quantum chaos system ||9|,ll0|] and its prediction has been 
recently confirmed by experiment. We derived an analytic optimal field to control the fully random matrix systems ITTll 
and the results are promising. However, the most realistic quantum systems may be modeled with a banded random 
matrix [5] and our previous method can not be directly applied to such a more general situation. 

In this contribution, we improve our previous approach 111 111 to deal with more realistic quantum systems with a 
banded random Hamiltonian. A new analytic field for optimal control is introduced and its validity is evaluated by 
numerically solving the Schrodinger equation for the multi-state-to-multi-state control problem. 

ANALYTIC EXTERNAL FIELD FOR CONTROLLING QUANTUM STATES 

In the previous works ITll[T2i[T3ll . we studied quantum control dynamics for a Hamiltonian driven by an external field 

e(0, 

H[e(t)]=H + e(t)V (1) 
where Ho and V are fully random matrices. For initial and target states in the eigenstate representation of Ho, 

I>i<p*) (2) 



((Pj\V\(p k ), Wjk = ^-J^ (3) 



l*oHE C >/>, |*r} = 
an analytic optimal field was derived as 
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where Vj k are elements for a fully random matrix with \V\ 2 = (|V;<fc| 2 ). Here () denotes an ensemble average. Driven 
by this field, the quantum state | y(t)) is shown to be steered from an initial state |<t>o) at t = to a target state |<J>r) at 
t = T according to 
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where Uq (t 2 , t\) is a propagator generated by Ho . 

We can summerize the reason why the analytic field works well for fully random matrix systems: 

• The random phase property of cj and dj is used, where the initial and target states are represented in linear 
combinations of many eigenstates \<Pj). If we introduce an approximate number of states, N ~ £ \cj\ ^ £ \d k \, 
sums of complex numbers, £c/ 2 , £of/ 2 , Y< c *jdj, etc., are quantities with an order 0(l/N), while £|c/| 2 = 
E|<//| 2 = 1 from the normalization condition. Thus, quantities of 0(1 /N) can be ignored compared to those 
of 0(1) forJV->°°. 

• The rotating wave approximation is applicable, which is valid when the control field amplitude is small enough. 
This situation is satisfied when T — > °°. 

Using these properties, we proved the validity of our analytic control field for fully random matrix systems |[Tll[T2i[T3ll . 



EXTENSION FOR BANDED RANDOM MATRIX SYSTEMS 

We consider the case that the interaction Hamiltonian V is a banded random matrix in the eigenstate representation of 
Hq. The elements of V are random complex numbers with distribution 
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We introduce an analytic optimal field as an extension of the analytic field for fully random matrix systems, 

e(t)=Y,^[Aj k c* j V jk d k e iw * t ], 
jk 

with an extra-amplitude factor A, The coefficients a At) satisfy the Schrtidinger equation 
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ih—a k (t) = [Aj k Cj*d k +A kj *c k dj*] \V jk \ aj(t) 
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under the rotating-wave approximation. If we assume that the transition is smooth, ai(t) should be written as 

a k (t) =c * cos (|f) - ^ sin (|f) ■ 
Substituting these coefficients into Eq. (0, we obtain a relation 
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under the assumption of random phases. Finally, we obtain conditions for the coefficients A i k 
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If we consider the case that those coefficients cj and d k have Gaussian distribution functions in the energy space, 
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FIGURE 1. Left: Schematic picture of the multi-state-to-multi-state control process. Right: An example of the interaction 
Hamiltonian V (128 x 128, Aq = 32), which is assumed to be a banded random matrix in the energy representation of Wq. 



with centers E c and E c / and widths A c and A4, we can define the analytic optimal field as 
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This field has a finite amplitude only when 



A c . < A and A d < A . 



(13) 



If not, the field has an infinite amplitude in the limit of Ej,E k — > ±°° by the exponential factor A: k in Eq. ( fT2l . Thus, 
the analytic field is refined when the widths of the initial and target states are relatively small compared to the width 
of the banded random matrix elements. 



Numerical evaluation 

We shall confirm the validity of the new optimal field for the system with a banded random-matrix interaction. The 
numerical test is configured as follows (see Figure[T|l. The initial and target states are defined as quantum vectors (O 
with random complex coefficients cj and di subject to ( fTTT i. Here, we choose A c = A c i = 32, E c = —16, and E c [ = 16, 
where Hq is a 128 x 128 random matrix of the Gaussian orthogonal ensemble and is scaled so that its eigenvalues {£; } 
are distributed in an interval [—64, 64]. The interaction Hamiltonian V is also a 128 x 128 matrix while its elements 
obey a banded-random distribution (0 in the eigenstate representation of Ho with Ao = 32. 

The optimal field ( flZb is calculated from those quantities {c,}, {d k }, {Vjk}, and {Ej} with parameters T and Ao. In 
order to check the validity of our optimal field (fT2l . we solve the initial value problems with Hamiltonian (fl]i driven 
by the optimal field ( fT2l for various band widths Ao of the interaction Hamiltonian V. The results are shown in Figure 
[2] When we use the original analytic field Q, the performance of the optimal field (dashed curve) decreases for the 
banded matrices with smaller widths. On the other hand, the final overlaps (solid curve) by the refined analytic field 
(fT2l does not change even for the smaller width untill the limit Ao w A c or Aq « A^. 




FIGURE 2. Numerical results for the multi-state-to-multi-state control problem. The solid (dashed) curve represents the result 
by the new optimal field Eq. J12I ) (by the previous optimal field with Ajf. = 1). (a) The final overlap is shown as a function of Ao 
(width of V). (b) The amplitude (defined as a time average of the absolute square of a field) of the optimal field as a function of Aq. 



CONCLUSION 

We extended our previous analytic approach for controlling complex quantum systems to deal with more realistic 
systems with a banded random matrix. The key ingredient is the amplitude factor A », which is an exponentially 
growing function, introduced in the analytic optimal field Eq. (fl2l i. We showed that the new analytic optimal field 
outperforms the previous optimal field for the multi-state-to-multi-state control problem. Interestingly we found a 
threshold of the width of the banded random Hamiltonian Ao for the control to be successful: Ao ~ A c , Aj where the 
latters are the width of the energy spreading of the initial and final states. In the near future we will apply this optimal 
field to quantum chaos systems such as quantum kicked rotors (tops) and to more realistic molecular systems 11411 . 
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